Magnetic- field dependence of dynamical vortex response in two-dimensional 
Josephson junction arrays and superconducting films 
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The dynamical vortex response of a two-dimensional array of the resistively shunted Josephson 
junctions in a perpendicular magnetic field is inferred from simulations. It is found that, as the mag- 
netic field is increased at a fixed temperature, the response crosses over from normal to anomalous, 
and that this crossover can be characterized by a single dimensionless parameter. It is described 
how this crossover should be reflected in measurements of the complex impedance for Josephson 
junction arrays and superconducting films. 

PACS numbers: 74.40. +k, 74.60. Ge, 74.76.-w, 75.40.Gb 
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Two-dimensional (2D) vortex physics is strongly re- 
flected in the properties of Josephson junction arrays 
(JJA) and supjerconducting films, as well as high- Tc su- 
perconductors .Ela In the absence of a magnetic field, the 
phase transition is of the Kosterlitz-Thouless (KT) type 
and is driven by the unbinding of thermally created 
vortex-antivortex pairs. In the presence of a perpendicu- 
lar magnetic field, ordered flux lattices are formed which 
melt at high enough temperatures. [Although these ther- 
modynamic aspects are well known,EI'n the dynamical re- 
sponse properties, which contain key information of the 
vortex physics and can be extracted from rneasurements 
of the complex impedance and the flux noise,EJ 13 are much 
less well understood.ElEI In the present paper we focus on 
the temperature region above the 2D flux lattice melting 
and infer that the dynamical response as a function of 
the perpendicular magnetic field has intriguing charac- 
teristics. 

The features of vortex dynamics are contained in the 
complex impedance Z{u)) or equivalently the conductiv- 
ity a{u!) = 1/Z{uj), which can be expressed as cr{uj) = 
l/i?Ar — l/iwLfee(a;) with the normal-state resistance Rn, 
the kinetic inductance L^, and the vortex dielectric func- 
tion l/e(a^|-|djescribing the effect of the vortices on the 
dynamics.tloO In the present work we obtain l/e{u;) and 
the resistance R — Re[a~^{uj — 0)] from simulations and 
characterize the vortex dynamics in terms of these quan- 
tities. 

We use the 2D resistively shunted junction (RSJ) 
model on an L X L triangular lattice with periodic bound- 
ary conditions (see Fig. m. The equations of motion in 
the presence of an external magnetic field B — Bz follow 
from the local current conservation 
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where Oi is the phase of the complex order parame- 
ter at site i, Gij is the lattice Green's function, and 
the primed summation is over six nearest neighbors of 
j. The thermal noise current rjjk in units of the sin- 
gle junction critical current Ic satisfies {rjijlt)) = and 
{'nij{t)'nki{0)) = 2T{S,k5ji - 6iiSjk)S{t), where (• • •) is the 



ensemble average and the temperature T is in units of 
the Josephson couphng strength J = UIc/2e. The mag- 
netic bond angle Aij = (27r/$o) /■"' A • dl with the flux 
quantum <I>o and the magnetic vector potential A = Bxy 
satisfies the plaquette sum ^ Aij — 27r/ with the frus- 
tration /, corresponding to the number of flux quanta per 
elementary triangle. The time t is in units of h/2eRNlc 
with the shunt resistance Rn, and in the following we 
simplify unit system to Rn = Ic = ?i/2e = 1. We use 
a second-order algorithm for numerical integration with 
the time step At — 0.05 and the system size L = 64. The 
key quantity is the time correlation function G{t), which 
is related with the dielectric function l/e(aj) by 
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and is determined from G{t) 
whej^po — \/3(cos(( 



(2poT/V3L2)(F(t)F(0)),^ 
Oj — Aij)) is the superfluid den- 



with 



sity,QH and F{t) = J2{ij) sin(6'i - 9j - Aij)e,j 
the sum over all links and ey being the unit vector from 
i to j. In our unit system po is the same as the kinetic 
inductance 1/Lfc so that the relation between cr(a;) and 
l/e{uj) reduces to 
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Figure pi shows the resistance i? as a function of frus- 
tration / at r = 0.1 and 0.5 obtained from G{t) using 
Eqs. (|) and (|) together with R = Re[l/a{uj = 0)]. In 
order to check the present method we have also calculated 
the resistance in a completely independent way using the 
fluctuating twist boundary conditions described in Rcf. 0. 
As seen in Fig. ^ both determinations give consistent re- 
sults but the present method appears to converge much 
better. The inset in Fig. g shows R versus T and verifies 
that T = 0.1 and 0.5 are both well above the resistive 
transition temperature. The naive expectation is that 
R is proportional to the [density of free vortices which in 
turn is proportional to /E This corresponds to a situation 
where all vortices effectively move diffusively as indepen- 
dent particles. As seen in Fig. this expectation is borne 



out for T = 0.5 for the two fowest frustratfons / = 1/64 
and 1/32 (these data points fall on a straight line through 
the origin). However, the highest frustration / = 1/16 
falls slightly below this line indicating a suppression of 
the free vortex diffusion. This effect is much more pro- 
nounced for the lower temperature T — 0.1, where both 
resistances for / = 1/16 and / = 1/32 fall well below the 
straight line through the origin. Since the-jiprtex pinning 
caused by a triangular array is very weak,E^ the essential 
part of the suppression is caused by some vortex correla- 
tion mechanism and not by a lattice pinning effect. 

Figure ^ shows the real and imaginary part of 1/€{lo) 
obtained from G(t) by aid of Eq. (0). One immediately 
notes a qualitative difference; in Figs, ^(a) and (b) the 
maximum of |Im[l/e(w)] is very close to the crossing 
point with Re[l/e(tj)], whereas in Fig. 0(c) this maximum 
is displaced to the right. In order to analyze these dif- 
ferences in a systematic way we parameterize the curves 
by an interpolation between two response forms. The 
first is the Drude form which corresponds to vortices ef- 
fectively diffusing as independent particles. This form is 
hence expected to apply to the region where the resis- 
tance is linear in / (as for / — 1/64 and 1/32 at T = 0.5 
in Fig. 0). The time correlation function in this case is 
of the exponential form G°™'^<'(i) = (2/OT)e-*'"o and the 
corresponding response function is given by 
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where e describes an effective static polarization cif the 
vortex interaction due to vortex-antivortex pairs Jj and 



Re[l/e(w)] 



for small 



The second response 



form describes the dynamical response of bound vortex- 
antivortex pairs whose phenomenological form has been 
proposed in Ref. |l|, Minnhagen phenomenology (MP), 
and is given by Ref. |: G^^{t) = (2/OT)[Ci(woO sinwoi- 
si{ijJot) cos uJot], where Ci(a;) = — J^ dtcost/t, si(x) = 
— J dtsint/t, and G^^ ex 1/t for large times. The cor- 
responding response function is given by 
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which goes to zero linearly with lu and gives a vanishing 
resistance. It is sometimes referred to as an anomalous 
vortex dynamics since it exhibits a nonanalytic divei:^ 
gence of the conductivity for small uj: (t{ljj) oc — Intj.tll 
This form has been shown, from simulations of various 
models of 2D vortex physics in the absence of magnetic 
field, ton idescribe the response of the low-temperature 
phasej3'0't3 where the dynamics is dominated by bound 
vortex-antivortex pairsE 

In general both bound pairs and free vortices can exist 
at the same time, e.g., for f — above the KT transition 
and for a finite / above the resistive transition, and the 
dynamical response is expected taryontain contributions 
from both of the response types .ErQ In the present work 



we use the interpolation introduced in Ref. ^i G'"*(i) = 
eG'^''"^°(t)G^^(t), which contains three parameters, e, 
(To, and luq, or equivalently e, r — ujo/ao, and R. Here r 
is the important parameter which changes the response 
from pure Drude for r = to pure MP for r — oo, 
and R satisfies ctq = R[pq{7t + rliir)]/{l — i?)[7re(l + r'^)] 
from G'"* together with Eqs. @ and (|). Since R can 
be independently determined without curve fitting as in 
Fig. ||, we are left with only two fitting parameters e and 
r. Furthermore, for a finite / above the resistive transi- 
tion there are few bound vortex-antivortex pairs present 
in the static limit so that e ~ 1,EJ and the dynamical 
response is essentially characterized by the single param- 
eter r. In the following analysis we use R in Fig. g and 
fit data to the dynamical response function l/e{uj) with 
two free parameters e and r. Figure shows three typical 
examples; the fits are very good and the obtained values 
e are in all cases close to 1 [ e w 1.0, 1.0, and 0.93 for 
Fig. 11(a), (b), and (c), respectively]. 

Figure H shows the characterization of the dynamical 
response in terms of the parameter r. The response for 
T = 0.5 corresponds to a small r and hence Drude-like 
behavior. This is consistent with independent diffusion 
of free vortices and with the linear relation between R 
and / found in Fig. 0. However, as seen in Fig. 0, r 
increases somewhat with /, which suggests an increase 
of "vortex-pair- like" correlations. This seems consistent 
with the small suppression of R with increasing / for 
T = 0.5 in Fig. 0. The same feature is much more dra- 
matic for r = 0.1: r increases from a small Drude-like 
value to a large MP-like value with /, which suggests that 
the response for the largest / is dominated by "vortex- 
pair-like" correlations, and is consistent with the large 
suppression of R at large / for T = 0.1 in Fig. g. 

An alternative way to analyze the crossover from 
Drude to MP response is to focus on the peak ratio, de- 
fined as the ratio |Re[l/e(tj)]/Im[l/e(ijj)]| at the maxi- 
mum of |Im[l/e(cj)]|. For the pure Drude response this 
ratio is unity whereas for the pure MP response it is 
2/7r « 0.64 [see Fig. g where (a) and (b) have peak ra- 
tios close to the Drude value whereas the peak ratio in 
(c) is close to the MP value] . Figure || displays the peak 
ratios obtained directly from the data (by obtaining a 
sufficiently number of data points close to the maximum) 
and shows that the crossover from Drude-like to MP-like 
behavior at T = 0.1 with increasing / is clearly signaled 
as a decrease in the peak ratio. 

Our simulations show that in a certain parameter 
range the dynamical response as a function of frustra- 
tion crosses over from a free-particle-like response to a 
response which has the same form as the response of 
bound vortex-antivortex pairs and that this crossover is 
linked to a suppression of the resistance. In Ref. the 
same type of crossover from Drude to MP response with 
increasing frustration was found for a 2D XY model with 
time-dependent Ginzburg-Landau (TDGL) dynamics. In 
Ref. O it was found that a triangular JJA with a small 
finite frustration shows a MP-like response well below the 



zero-frustration KT transition and in Ref. it was sug- 
gested, based on the peak ratio, that the same is true for 
the superconducting film measured in Ref. ^ 

Several theoretical attempts have been ma dij , t o under- 
stand the origin of the anomalous dynamics :tillla Ref. 04 
ascribes the effect to a coupling between spin waves and 
vortices, suggesting that there should be anomalous dy- 
namics for TDGL but not for RSJ, in contrast to what 
has been demonstrated in the present work. Ref. uB at- 
tributes the effect to vortex correlations which should 
apply to both the TDGL and RSJ. However, none of 
the theoretical formulations yields a crossover to the MP 
form and a shift of the peak ratio from 1 to 2/7r with 
increasing magnetic field and hence with an increasing 
resistance. Thus an adequate theoretical explanation is 
still largely lacking. In Ref. it was found from simula- 
tions that an increase of the magnetic field also causes an 
increased density of antivortices suggesting that vortex- 
antivortex correlations may play a role with a possible 
connection to the fact that the MP response also de- 
scribes the response of the low-temperature phase for 
/ = 0. 

Our analyzing method makes it possible to assess to 
what extent the crossover from normal to anomalous re- 
sponse described here shows up in measurements on real 
systems. 
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FIG. 2. Resistance _R as a function of frustration / at 
T — 0.5 (squares) and 0.1 (circles). The empty symbols rep- 
resent R obtained as described in text, and the filled symbols 
as obtained by the method described in Ref. q (denoted by 
Ra)- The dashed lines through the origin correspond to ex- 
pectation from free vortex diffusion. The full drawn lines are 
guides to the eyes. The inset shows that the onset of resis- 
tance starts well below T = 0.1. 
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FIG. 4. The parameter r characterizing the response as a 
function of /. At T = 0.5 r increases shghtly with increasing / 
but remains much smaUer than 1, characterizing a Drude-hke 
response. At T = 0.1, r changes from r ^ 1 to r 3> 1 with in- 
creasing /, signifying a crossover from Drude-hke to MP-Iike 
response. 
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FIG. 3. The response function l/e(cij) for (a) T = 0.5 and 
/ = 1/16, (b) T = 0.1 and / = 1/64, and (c) T = 0.1 
and / = 1/16. The full and dashed lines are the fits to the 
functional form described in text. For each pair (T, /) this 
essentially gives a one-parameter characterization. 




FIG. 5. Peak ratio as a function of frustration /. At 
T = 0.5 the peak ratio decreases slightly with increasing / 
but remains close to the Drude value 1. At T = 0.1 the peak 
ratio crosses over from 1 to the MP value 2/7r. 



